基于 交叉 点 数据 和 三 维 非 参数 模型 的 雷达 高 度 计 海 况 偏差 估计 方法 
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摘 要 : bres s State Bias, SSB) 是 雷达 高 度 计 测量 海面 高 度 的 重要 误差 源 。 目 前 ， 业 务 化 运行 的 雷达 高 度 计 
的 海 况 偏差 校正 都 是 采用 基于 风速 (中 和 有 效 波 高 (SWH) 的 2 维 经 验 模型 方法 , 其 海 况 偏差 校正 不 确定 度 约 为 2 cm. 
该 文 提 出 了 一 种 基于 交叉 点 数据 和 3 维 非 参 数 模型 的 海 况 偏差 估计 方法 ， 该 方法 利用 星 下 交叉 点 数据 ， 采 用 基于 
U, SWH 以 及 平均 波 周 期 (MWP) 的 3 维 非 参 数 模 型 进行 海 况 偏差 估计 。 该 文 利 用 这 种 估计 方法 对 Jason-2 卫星 雷达 
高 度 计 2009~2011 三 年 的 数据 进行 了 处 理 ， 将 处 理 结果 与 Jason-2 卫星 高 度 计 的 地 球 物 "1 3 8 E (Geophysical Data 
Records, GDR) 中 的 海 况 偏差 校正 项 进行 对 比 ， 结 果 表 明 该 文 提出 的 海 况 偏差 估计 方法 平均 能 降低 1.64 cm? 的 交叉 
点 海面 高 度 不 符 值 的 方差 和 0.92 cm? 的 沿 轨 海 面 高 度 异常 值 的 方差 ,分 别 对 应 于 1.28 cm 和 0.96 cm 的 均 方 根 (RMS) 
海面 高 度 ， 这 对 于 提高 高 度 计 数据 产品 的 精度 具有 重要 意义 。 
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Sa Abstract: The Sea State Bias (SSB) is an important source of error in satellite altimetry. Operational SSB correction 


models are based on the altimeter-measured Wind Speed (U) and Significant Wave Height (SWH). This paper presents a 


new method to estimate the SSB from the crossover differences using a three-dimensional nonparametric model based on 
U, SWH and the Mean Wave Period (MWP). Evaluated by the separate annual data sets from 2009 to 2011, the SSB 
values estimated with the presented method can decrease the variance of the crossover Sea Surface Height (SSH) 
differences by 1.64 cm?, or 1.28 cm RMS, and decrease the variance of the Sea Level Anomalies (SLA) by 0.92 cm? , or 
0.96 cm RMS in comparison to the SSB values in the Geophysical Data Records (GDR) of Jason-2. It is of great 
significance for improving the precision of altimeter products. 
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全 球 海平 面 上 升 已经 越 来 越 引 起 人 们 的 关注 ， 雷 达 高 度 计 的 一 个 主要 应 用 就 是 测量 平均 海面 高 度 [1-3]。 
i 摆 是 非 高 斯 分 布 的 ， 波 谷 的 曲率 半径 比 波峰 大 ， 所 以 波 谷 能 反射 更 多 的 电磁 能 量 ， 从 而 使 得 
度 计 测 得 的 平均 海面 低 于 真实 的 平均 海 ， 这 种 效应 称 为 电磁 偏差 中。 在 高 度 计 的 回 波 跟 踪 算 法 中 ， 通 
段 设 海面 高 度 的 概率 密度 函数 为 高 斯 函数 ， 而 实际 海面 高 度 的 概率 密度 函数 是 有 偏 斜 度 的 ， 这 样 就 会 产 
扁 斜 度 偏 差 。 电 磁 偏 差 和 偏 斜 度 偏差 一 起 统称 为 海 况 偏差 (Sea State Bias, SSB), 海 况 偏差 的 不 确定 度 很 大 ， 
够 达到 2 cmrzs o， 随 着 精密 定 轨 技 术 的 发 展 ， 在 Jason 系列 测 高 卫星 中 ， 海 况 偏差 已 经 取代 轨道 误差 成 为 
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圭 达 高 度 计 测 高 最 大 的 误差 源 [6 站 。 
国际 上 主要 从 理论 模型 G1 和 经 验 模 型 [5 两 个 方面 对 海 况 偏差 进行 研究 。 海 况 偏 差 的 机 理 非常 复 
杂 ， 目 前 理论 模型 并 不 实用 ， 业 务 化 运行 的 高 度 计 都 是 采用 经 验 模 型 对 海 况 偏差 进行 校正 。 海 况 偏差 经 验 
模型 包括 参数 模型 [i631 和 非 参数 模型 [512.1410J， 非 参数 模型 精度 更 高 ， 被 广泛 应 用 于 海 况 偏差 校正 。 传 统 
的 海 况 偏差 经 验 模 型 使 用 高 度 计 测 得 的 U 和 SWH 作为 模型 的 输入 ， 但 仅仅 由 U RI SWH 只 能 得 到 有 限 的 
海 况 偏差 信息 , 新 参量 的 引入 成 为 海 况 偏 差 经 验 模 型 新 的 研究 方向 G5160]。 文献 [15,16] 在 2010 年 利用 全 球 波 
模式 数据 WAVEWATCH III 生成 了 平均 波 周期 (MWP), 利用 海面 高 度 异常 数据 建立 了 一 个 基于 U, SWH 和 
MWP 的 3 维 海 况 偏差 估计 模型 ， 结 果 表 明和 传统 的 基于 U 和 SWH 的 e 维 模 型 相 比 ，s 维 海 况 偏差 估计 模 
型 可 以 减 小 1.26 cm? Bii ie BEJI 25002. AA 3 维 海 况 模 型 中 的 平均 海平 面 是 通过 不 同 高 度 计 多 年 测 高 数 
据 平均 得 到 ， 这 样 会 引入 较 大 误差 ， 由 此 得 到 的 SSB 估计 值 误 差 也 比较 大 。 
在 文献 [15,16] 的 基础 上 ， 本 文 提 出 了 一 种 基于 交叉 点 数据 和 3 维 非 参数 模型 的 海 况 偏差 估计 方法 。 本 
文 利 用 非 线性 海浪 理论 简单 分 析 了 波 周 期 和 海 况 偏差 的 大 致 关 系 ， 为 将 波 周期 信息 引入 海 况 偏差 估计 提供 
了 理论 依据 。 利 用 本 文 方法 分 别 对 Jason-2 卫星 雷达 高 度 计 2009-2011 三 年 的 数据 进行 了 处 理 ， 并 将 处 理 结 
果 与 Jason-2 卫星 高 度 计 的 GDR 产品 中 的 海 况 偏差 校正 项 进行 对 比 ， 结 果 表 明 本 文 提出 的 海 况 偏差 估计 方 
法 能 极 大 地 降低 交叉 点 海面 高 度 不 符 值 的 方差 和 沿 轨 海面 高 度 异 常 值 的 方差 ， 这 对 于 提高 雷达 高 度 计 的 测 
高 精度 具有 重要 意义 。 
2 波 周期 与 海 况 偏差 的 关系 

文献 [15,16] 曾 将 平均 波 周期 引入 海 况 偏差 估计 中 ， 但 是 国内 外 还 没有 公开 资料 论述 波 周期 与 海 况 偏差 
的 关系 ， 本 文 希 望 通过 海浪 理论 将 波 周 期 与 海 况 偏差 联系 起 来 ， 从 而 为 将 波 周期 信息 引入 海 况 偏差 估计 提 
供 理论 依据 。 海 况 偏差 主要 是 由 海面 的 非 线性 特性 所 引起 的 ， 所 以 可 以 使 用 非 线性 海浪 理论 来 进行 研究 。 
常用 的 非 线 性 海浪 理论 主要 有 斯 托 克 斯 (Stokes) 波 理论 、 椭圆 余弦 (Cnoidal) 波 理论 和 孤 江 (Solitary) 波 理论 , A6 
圆 余 弦 波 理论 、 扳 立波 理论 以 及 五 阶 Stokes 波 理论 都 主要 适用 于 浅水 ， 而 高 度 计 的 绝 大 部 分 有 效 测量 点 都 
在 深水 区 ， 所 以 本 文 使 用 四 阶 Stokes 波 理论 。 
根据 四 阶 Stokes 波 理论 ， 波 面 轮廓 可 以 表示 为 cl5]: 
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记分 别 为 波长 、 波 高 和 波 数 ，T 为 波 周期 ，< 为 与 波 
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AF, L= Z a+ ka’), H= 2a(l+ Ska’) 和 k = 


p L 
高 及 有 关 的 待定 实 常 数 。 波 面 轮 廊 V (x) 满足 : 
à, VG- Ga? + kia) = 0 (2) 


RORI, z, = Ska? + Kia 为 静止 时 的 自由 面 ， 或 者 说 坐标 原点 位 于 静止 水 面 以 下 3ka? + k'a* 处 ， 
波峰 到 坐标 原点 的 距离 为 


H = Dr + ka + a(0 3 ka?) + ir + i8 
^ 2 8 2 6 


+ “ka? 十 二 KG 
8 3 
波 谷 到 坐标 原点 的 距离 为 : 
H = iu + k/a^- a(14 2 kg?) 4 (ip 十 i55 
"73 8 2 6 (4) 


" pr + laa 
8 3 


半 波 面 高 度 到 原点 的 距离 为 : 
H, = liu + H,)= Lg? fat (La? 十 pego. Ip (5) 
h 2 c 2 2 6 3 
半 波 面 高 度 到 静止 水 面 的 距离 为 : 


1 11 1 


H, = H, = Zo = P + r3 十 rid (6) 
式 (6) 表 明 水 面 质 点 的 振动 中 心 不 在 静 止 水 面 上 ， 而 是 高 于 静止 水 面 ， 高 出 静止 水 面 的 距离 约 为 H,， 


zo = Ska? + ki 可 以 看 作 是 高 度 计 测 得 的 平均 海面 ， 太 , 可 以 看 作 是 真实 的 平均 海面 ， 也 就 是 高 度 计 测 得 


的 平均 海面 会 低 于 实际 的 平均 海面 ， 而 真实 的 平均 海面 与 高 度 计 测 得 的 平均 海面 之 差 就 可 以 看 作 是 电磁 偏 
差 ， 所 以 有 H, 也 可 以 看 作 是 四 阶 Stokes 波 理论 得 到 的 电磁 偏差 。 
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e 图 1. 四 阶 Stokes 波 理论 得 到 的 电磁 偏差 H, Bic H A 并 的 变化 
图 1(a) 给 出 了 当 波 周期 工 固定 为 9 s IM. H 随 波 高 五 的 变化 关系 ， 当 波 周期 T JEDER. H, 随 波 高 H 
C5) ”的 增 大 而 增 大 ; 图 1(b) 给 出 了 当 波 高 H EEA smit, H 随 波 周期 的 变化 关系 ， 当 波 高 HEN, H, 
N 随 波 周期 T AKTRA. M ERE ERARE RRRS MAR 1 也 可 以 大 致 反映 海 况 偏差 与 波 
: 长 和 波 周 期 的 关系 。 
以 上 提 到 的 波浪 是 具有 固定 波 高 、 周 期 的 单一 波浪 ， 其 波 面 形状 是 规则 的 ， 然 而 实际 海面 上 的 波浪 都 


是 不 规则 的 ， 由 不 同 波 高 和 周期 的 单一 波浪 舍 加 而 成 。 所 以 ， 在 利用 经 验 模型 来 佑 计 海 沈 偏 差 时 ， 我 们 使 
用 有 效 波 高 和 平均 波 周 期 。 
3 已 有 的 s 维 非 参 数 海 况 偏差 模型 
3.1 3 维 非 参数 回归 模型 
假设 > 为 因 变量 ，x = (x,x,,x,) 为 自 变 量 ， 那 么 一 般 的 3 维 回 归 模 型 可 以 表示 为 : 

y= r(x)+e (7) 
AP ra) AEIR e 为 零 均值 的 误差 项 ， 即 E(e)= 0 。 给 定 n 组 观测 值 Q ,x,),i = L2..sn, Ex 的 邻 
域 对 回归 水 数 r(x) 进行 一 阶 Taylor EFAA: 


rz) = bt & bZ- x) (8) 


式 中 ，z 是 一 个 在 x 的 邻 域 含有 s 个 变量 的 向 量 。 通 过 找到 一 组 = (usus) 使 得 式 (9) 最 小 : 


Bb,- po- beg x fce x) , 
这 是 典型 的 加 权 最 小 二 乘 回归 问题 ，b 的 估计 值 为 
B= XIWX,) 'XIWy ka 
ui Xa X, Xu 7X, 
AF, X, = : e m , W = diag{K s(x -x,)} 是 一 个 (n”n) 的 加 权 和 矩阵 ，K slx -x,) 通 


常 是 由 核 函 数 构成 的 概率 密度 函数 。 

r(x) 的 局 部 线性 Fi au Go) Ju, , TUB 的 其 它 成 分 是 函数 r 的 一 阶 偏 导数 ， 即 

Pus) $ o eT QCIWX y X Wy (11) 

AF, y'= D ,yy，]， e 是 一 个 单位 可 量 ， ei = [1,0,0,0]. 
3.2 文献 [15,16] 的 3 维 非 参数 海 况 偏差 模型 

传统 的 非 参 数 海 况 偏差 模型 主要 是 利用 U f SWH 作为 模型 输入 的 o 维 经 验 模型 ， 文 献 [15,16] 在 2010 
年 提出 的 s 维 非 参数 海 况 偏差 模型 采用 U, SWH 和 MWP 作为 模型 的 输入 De。 如 何 提取 SSB 信息 是 建立 SSB 
模型 的 关键 ，TRAN 等 人 将 未 经 SSB 校正 的 瞬时 海面 高 度 减 去 通过 平均 海面 模型 得 到 的 平均 海面 高 度 ， 从 
而 得 到 未 经 SSB 校正 的 海面 高 度 异常 值 并 从 中 提取 SSB 信息 。 

假设 海 况 偏差 SSB 是 x = (U,SWH, MWP) 的 函数 g ， 即 SSB = gx), SSH, 表示 未 经 SSB 校正 的 海面 
高 度 , MSS 表示 由 平均 海面 模型 得 到 的 平均 海面 , y, = SSH, - MSS, 表示 未 经 SSB 校正 的 海面 高 度 异常 值 ， 
由 文献 [15] 和 文献 [16] 可 知 ， 观 测 方程 可 以 写 为 ， 
y, = (SSB), +e, = 8Cc)+ei (12) 
这 样 就 得 到 了 式 (7) 中 典型 的 回归 模型 的 表达 式 ， 给 定 n 组 观测 值 Q,,x,) ， 通 过 式 (11) 就 可 以 得 到 任意 xz 所 对 
应 的 SSB fie Go) 。 
4 基于 交叉 点 数据 的 3 维 非 参 数 海 况 偏差 估计 
4.1 基于 交叉 点 数据 的 3 维 非 参数 海 况 偏差 估计 方法 

与 文献 [15,16] 提 取 SSB 信息 的 方法 不 同 , 本 文 从 未 经 SSB 校正 的 交叉 点 海面 高 度 不 符 值 中 提取 SSB f 
E. RERUM RO 

SSH 表示 未 经 SSB 校正 的 海面 高 度 , 在 交 义 点 处 , 通过 降 轨 和 升 轨 的 SSH 相 减 ， 可 以 消除 大 地 水 准 面 
中 不 随时 间 变 化 的 部 分 ， 假 设 SSB 是 x = (U,SWH,MWP) 的 函数 ge ， 即 SSB = g(x)， 令 y = SSH,-SSH,; 
根据 文献 [5] 和 文献 [10]， 观 测 方程 可 以 写 为 : 


c 


y= gx)- gx)*e (18) 


在 给 定 x, =x 的 条 件 下 ，y 的 条 件 期 望 为 : 
Ely |x, 2 x]2 g(x)- Elgx)|x, = x] (14) 
根据 非 参 数 估计 理论 ， 条 件 期 望 rc)= Ely |x] 的 核 估计 器 Bx) 可 以 看 作 是 一 系列 观测 值 y 的 加 权 平 均 : 


$a) - a yax- x,) (15) 


BUR La, (x - x,) 可 以 通过 式 (11) 得 到 。 在 式 (11) 中 ， 令 T = e A WX, XW ， 可 以 得 到 权 函 数 : 
a (x- x,)=T() (16) 
给 定 n 组 观测 值 0,,x,,x,) ， 利 用 式 (14) 和 式 (15) 可 以 得 到 : 


ge)= Nba Ge xe qeu G x) (17) 
只 要 知道 了 g(x,) ， 就 可 以 得 到 任意 x 所 对 应 的 gr) ， 将 式 (17) 中 的 x 蔡 换 为 r，， 可 以 得 到 


gx,)- Jia, Qr, ` xy) + gx a, CE " X) (18) 
i pcd. n 
转化 为 矩阵 的 形式 可 以 得 到 : 
(I -A)g = Ay (19) 


NPI Æa’ n) WAER, A 9E n n DIEE, CRA aua ux) e 


8; = [gE 8E a) 8E) y = yy yl 


8 ,不 能 通过 线性 系统 式 (19) 直 接 求 得 ， 因 为 7 -4 是 一 个 奇异 矩阵 ， 它 的 秩 为 n-1， 为 了 消除 这 种 不 确定 
性 ， cce. 个 元 素 强 加 一 个 值 ，g(x,,) = g。， 式 (19) 可 以 重 写 为 : 
Bg-Ay- Bg, (20) 
式 中 g 包含 g | 中 除 gC ) ZA n-i 427638. g^ = [g(x,,),…,8(x,,)]，B, 是 I -A 的 第 1 列 ，B, 是 其 余 列 ， 
g 的 最 小 二 乘 解 为 : 


(BTB ) BL (Ay - B,g) 21) 


再 联合 g(x,) = g,， 就 可 以 得 到 g, ， 将 它 带 入 式 (17) 就 可 以 得 到 任意 x 所 对 应 的 g(x) Bh Bo) 。 
对 于 每 一 个 输入 变量 x ， 要 得 到 其 对 应 的 SSB 估计 值 都 需要 经 过 大 量 的 矩阵 运算 ， 这 很 难 直 接应 用 于 
高 度 计 数据 的 实时 处 理 当 中 ， 为 了 解决 这 个 问题 ， 可 以 制作 一 个 SSB 估计 值 关 于 U, SWH 和 MWP 的 3 维 
查找 表 ， 分 别 计 算 查 找 表 中 每 个 网 格 点 的 输入 变量 x 所 对 应 的 SSB 估计 值 ， 再 通过 插值 的 方法 得 到 任意 高 
度 计 测量 点 的 SSB 估计 值 。 
4.9 两 种 3 维 非 参数 海 况 偏差 估计 方法 的 比较 
利用 海面 高 度 异 常数 据 估 计 SSB 和 利用 交叉 点 数据 估计 SSB 的 主要 区 别 是 获取 SSB 信息 的 方式 不 同 ， 
前 者 是 从 海面 高 度 异 常 值 中 获取 SSB 信息 ， 而 后 者 是 从 交叉 点 的 海面 高 度 不 符 值 中 获取 SSB 信息 。 与 利用 
海面 高 度 异 常数 据 估计 SSB 的 方法 相 比 ， 利 交叉 点 数据 估计 SSB 的 方法 具有 以 下 优势 : 
(1) 海 面 高 度 异常 值 是 高 度 计 测 得 的 瞬时 海面 高 度 与 通过 平均 海面 模型 计算 得 到 的 平均 海面 高 度 之 差 ， 
对 于 同一 个 平均 海面 模型 ， 在 某 一 个 位 置 ， 平 均 海面 高 度 是 固定 的 ， 其 中 高 度 计 的 系统 误差 会 影响 SSB 的 
测量 结果 ， 由 于 高 度 计 的 系统 误差 在 交叉 点 数据 获取 过 程 中 是 不 变 的 ， 通 过 在 交 又 点 处 将 升 轨 和 降 轨 的 海 
面 高 度 做 差 就 可 以 消除 这 样 的 误差 。 
(2) 计 算 海 面 高 度 异常 值 需 要 用 到 平均 海面 模型 ， 平 均 海面 模型 是 通过 多 个 高 度 计 多 年 测 高 数据 平均 得 
». 到 ， 不 同 高 度 计 所 采用 的 数据 处 理 方 法 有 差异 ， 所 用 到 的 误差 校正 方法 也 不 一 样 ， 这 些 都 会 引入 误差 。 同 
c 时 ， 利 用 不 同 平均 海面 模型 得 到 的 SSB iE ARAK, 而 利用 交叉 点 的 方法 则 没有 这 一 问题 。 
à 5 Jason? 卫星 高 度 计 海 况 偏差 估计 
e 5.1 交叉 点 数据 处 理 
CN 本 文 使 用 的 平均 波 周 期 来 自 于 欧洲 中 期 天 气 预 报 中 ， 人 Centre for Medium-range Weather 
~y Forecasts, ECMWTF) 再 分 析 数 据 集 ERA-Interim, ERA-Interim 是 ECMWF 的 第 s 代 全 球 再 分 析 数 据 产品 ， 
vc» RETE 1979 年 以 来 的 再 分 析 资 料 ， 并 实时 更 新 。ERA-Interim 提供 每 天 4 个 时 刻 的 数据 ， 分别 在 00: 00, 
^em 06: 00, 12:00 和 18:00 UTC， 数 据 以 网 格 的 形式 提供 ， 本 文 使 用 数据 的 网 格 分 辨 率 为 0.25 按 0.22. 
ERA-Interim 只 提供 每 天 4 个 时 刻 的 数据 ， 所 以 需要 将 数据 在 空间 和 时 间 上 插值 到 高 度 计 测量 点 处 。 本 文采 
一 用 先 在 空间 上 双 线 性 插值 ， 再 在 时 间 上 线性 插值 的 方法 。 本 文 用 到 的 风速 、 有 效 ; M 以 及 计算 海面 高 度 所 
FI 需要 的 数据 都 来 源 于 Jason-2 卫星 经 过 完全 定 标 的 D 版 本 的 GDR 数据 。Jason-2 卫星 于 2008 年 6 月 发 射 ， 
它 是 TOPEX/POSEIDON, Jason-1 测 高 卫星 的 后 继 星 ， 被 公认 为 目 TARRAT 的 卫星 。 
利用 交叉 点 数据 和 3 维 非 参 数 模型 估计 SSB 时 需要 用 到 风速 、 有 效 波 高 、 平 均 波 周期 和 交叉 点 处 未 经 
SSB 校正 的 海面 高 度 不 符 值 ， 需 要 把 它们 插值 到 交叉 点 处 。 本 文采 用 三 次 样 条 插值 的 方法 ， 只 有 在 交叉 点 两 
侧 各 有 4 个 连续 的 测量 点 时 得 到 的 交叉 点 才 有 效 ， 这 样 可 以 保证 较 好 的 高 度 计 数据 质量 
5.2 Jason-2 海 况 偏差 估计 过 程 
在 利用 非 参数 估计 方法 估计 SSB 时 ， 需 要 建立 一 个 关于 U, SWH 和 MWP 的 3 维 查 找 表 ， 图 2 给 出 了 
U, SWH 和 MWP 的 分 布 直方 图 , 图 中 的 结果 是 利用 2011 年 全 年 的 数据 统计 得 到 的 。 在 本 文 所 采用 的 查找 
表 中 ，U, SWH 和 MWP 的 取 值 范围 分 别 为 : 0-30 m/s, 0~12 m 和 0~18s。 在 利用 2 维 非 参数 模型 估计 SSB 


时 ， 在 得 到 B(x) 后 ， 需 要 整体 减 去 8 = 0SWH = 0, 4U = 0,SWH = 0) = 0 ， 即 平坦 无 风 的 海面 的 
海 况 偏差 为 零 。 在 利用 s 维 非 参数 模型 估计 SSB H, EAA, KLARAH io 整体 减 去 
$u = 0, SWH = GMWP = 9s) 的 方法 得 到 最 终 的 SSB 估计 值 的 查找 表 , 选择 MWP=9 s 是 因为 MWP 的 平 
均值 在 o s 附近 。 利 用 交叉 点 数据 和 维 非 参数 模型 估计 SSB 的 步骤 如 下 ; 
步骤 1 对 单个 重复 周期 (Cycle) 的 交叉 点 数据 ， 假 设 一 个 gGr，) go g 可 以 是 任意 一 个 合理 的 值 ， 
例如 sg， = - 10cm。 
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步骤 o 通过 式 (21) 得 到 WA g, 得 到 8 . 


步骤 3 通过 式 (17) 计 算 g(z) 在 (U，SWH，MWP) 的 三 维 网 格 的 值 ， 得 到 一 个 3 维 查找 表 。 
步骤 4 对 每 个 Cycle 的 交叉 点 数据 重复 步骤 1， 步 又 2， 步 又 3， 再 对 所 有 Cycle 得 到 的 结果 取 平 均 得 


到 g(x)。 


33k 5 在 所 有 网 格 点 , 将 计算 得 到 的 g(x) 同时 减 去 8C = 0SWH = 0,MWP = 95) 就 得 到 s 维 查找 表 


中 每 个 网 格 点 所 对 应 的 SSB 估计 值 。 


步骤 6 在 得 到 3 维 SSB 查找 表 后 ， 根 据 沿 轨 任 意 一 个 测量 点 的 U，SWH 和 MWP， 通 过 插值 的 方法 


得 到 该 测量 点 的 SSB 估计 值 。 


mean = 7.73 m/s 
std = 3.51 m/s 
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5.3 Jason-2 海 况 偏差 估计 结 


分 别 使 用 了 2009, 2010 和 2011 三 年 的 交叉 点 数据 和 3 维 非 参 数 模型 得 到 SSB 估计 值 (CR3D), 因为 利用 
不 同年 份 数据 得 到 的 SSB 估计 值 相差 不 大 ， 在 数据 密集 区 都 在 1 cm 以 下 ， 这 里 仅 以 利用 2011 年 的 数据 得 
到 的 SSB 估计 值 的 查找 表 作 为 参考 。 为 了 更 好 地 表示 结果 ， 固 定 第 3 个 输入 参量 ， 以 2 维 网 格 的 形式 给 出 
SSB 估计 值 在 各 个 平面 的 等 值 线 分 布 , 如 图 3 所 示 , 等 值 


的 平均 值 附 近 ， 这 样 可 以 保证 在 相应 的 平面 有 较 多 的 


图 2。 风速 U、 有 效 波 高 SWH 和 平均 波 周 期 MWP 的 分 布 直 方 图 


5 10 15 20 
MWP (s) 


线 的 单位 为 cm, 3 个 固定 值 分 别 选 在 U, SWH, MWP 
测量 点 ，3 个 固定 值 分 别 为 : MWP=9 s, Uz7.5m/s, 


SWH-2.5 m。 图 中 阴影 部 分 为 测量 点 个 数 小 于 20 的 区 域 ， 代 表 数 据 稀疏 区 域 ， 
集 区 域 ， 统 计数 据 来 源 于 2011 年 全 年 的 数据 ， 我 们 主要 关心 的 是 数据 密集 区 域 。 

图 3(a) 是 在 MWP 等 于 9s 时 , SSB 估计 值 在 (U, SWH) 平 面 的 分 布 图 , 当 MWP 固定 时 , 和 器 相 比 , SSB 
随 SWH 的 变化 更 为 明显 ， 在 U 给 定 的 条 件 下 ，SSB 的 幅 值 是 SWH 的 递增 函数 ， 这 和 传统 的 利用 基于 U 
和 SWH 的 2 维 非 参 数 模型 得 到 的 结果 5512, 二 是 一 致 的 ， 也 和 图 1(a) 的 结果 吻合 。 图 3(b) 是 当 U SET. 7.5 m/s 
ET, SSB 估计 值 在 (MWP SWH) 平 面 的 分 布 图 ， 当 UF SWH 固定 时 ，SSB 会 随 着 MWP 变化 ,，MWP RAK, 
SSB 的 幅 值 越 小 。 高 度 计 的 绝 大 部 分 有 效 测量 点 都 在 深水 区 ， 而 在 深水 区 ， 影 
在 波 高 一 定时 ， 波 周期 越 大 ， 波 陡 越 小 ， 海 面 非 线性 越 弱 ， 海 况 偏差 也 就 越 小 


也 和 图 1(b) 的 结果 吻合 。 在 数据 密集 区 域 ， 对 于 给 定 


的 SWH， 随 着 MWP 的 


而 图 中 白色 部 分 代表 数据 密 


响 海面 非 线 性 的 主要 是 波 陡 ， 
， 所 以 得 到 的 结果 是 合理 的 ， 
变化 ，SSB 有 约 3~4 cm 的 变 


化 。 图 3(c) 给 出 了 当 SWH 等 于 2.5m 时 ，SSB 估计 值 在 (U, MWP) 平 面 的 分 布 图 , 在 数据 密集 区 域 , 当 SWH 
固定 时 ，SSB 随 MWP 的 变化 比 随口 的 变化 更 为 明显 ， 这 表明 MWP 对 SSB 的 影响 比 上 对 SSB 的 影响 大 ， 


在 其 它 条 件 都 相同 时 ，MWP 越 大 ，SSB 的 幅 值 越 小 。 


MWP (s) 


(a) MWP —9s (b) U = 7.5 m/s 


4 6 8 10 12 14 16 18 
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(c) SWH = 2.5 m 


图 3 利用 CRsD 得 到 的 Jason-2 SSB 估计 值 (cm) 在 各 平面 的 分 布 


表 1 利用 交叉 点 数据 和 3 维 非 参 数 模型 得 到 的 SSB 估计 值 与 Jason-2 的 GDR 中 的 SSB 校正 项 的 比较 


ASSH 方差 (cm?) SLA 方差 (cm?) 
SSB fii H8 验证 数据 集 验证 数据 集 

2009 2010 2011 2009 2010 2011 
JA2GDR 25.16 25.44 25.39 44.66 46.21 49.64 
2009 version, CR3D 23.629 23.838 23.85 43.95 45.39 48.45 
2009 version, difference (CR3D - JA2GDR) -1.54 -1.61 -1.54 -0.71 -0.82 -0.79 
2010 version, CR3D 23.53 23.71 23.14 43.79 45.13 48.54 
2010 version, difference (CR3D — JA2GDR) -1.63 -1.73 -1.65 -0.94 -1.08 -1.10 
2011 version, CR3D 23.52 23.69 23.71 43.83 45.26 48.50 
2011 version, difference (CR3D — JA2GDR) -1.64 -1.75 -1.68 -0.83 -0.95 -1.04 


注 : 表 中 ASSH 方差 表示 交叉 点 海面 高 度 不 符 值 的 方差 ,第 1 栏 中 的 JA2GDR” 表 示 Jason-2 的 GDR 产品 中 的 SSB 校正 项 ,“2009 version, CR3D” 
表示 利用 2009 年 的 交叉 点 数据 和 3 维 非 参数 模型 得 到 的 SSB 估计 值 , “2009 version, difference (CR3D - JA2GDR)" 所 在 的 行 表 示 相 对 于 Jason-2 
的 GDR 产品 中 SSB 校正 项 ， 利 用 2009 年 的 交叉 点 数据 和 3 维 非 参数 横 型 得 到 的 SSB 估计 值 所 计算 得 到 的 ASSH 或 SLA 的 方差 的 减少 量 。 


Jason-2 的 GDR 产品 中 含有 SSB 校正 项 ， 本 文 将 利用 交叉 点 数据 和 3 维 非 参数 模型 得 到 的 SSB 估计 值 
(CR3D) 和 Jason-2 的 GDR 产品 中 SSB 校正 项 (JA2GDR) 作 比较 ， 以 此 来 评估 利用 本 文 所 提出 的 SSB 估计 方 
法 得 到 的 SSB 估计 值 。 由 于 目前 还 不 能 直接 利用 卫星 高 度 计 自身 测 得 的 数据 得 到 海 况 偏差 的 精度 ， 所 以 通 
常 是 通过 海 况 偏差 校正 对 整个 测 高 系统 的 贡献 来 评估 SSB 估计 值 。 

海面 高 度 不 符 值 是 最 主要 的 用 于 分 析 整 个 测 高 系统 性 能 的 工具 ， 通 过 海面 高 度 不 符 值 可 以 分 析 升 轨 和 
降 轨 的 海面 高 度 (Sea Surface Height, SSH) 在 交叉 点 处 的 一 致 性 。 在 理想 的 条 件 下 ， 交 叉 点 海面 高 度 不 符 值 应 
该 为 零 ， 但 是 由 于 轨道 误差 、 地 球 物理 校正 误差 和 海洋 波动 等 因素 的 影响 ， 交 叉 点 海面 高 度 不 符 值 通常 不 
wi 为 零 。 交 叉 点 海面 高 度 不 符 值 的 方差 可 以 估计 整个 测 高 系统 的 表现 ， 其 值 越 小 ， 表 明 测 高 系统 的 整体 表现 
n BDs, MAR 1 可 以 看 出 ， 和 Jason-2 的 GDR 中 的 SSB 校正 项 相 比 ， 在 其 它 条 件 都 相同 时 ，CRsD 可 以 
C 将 海面 高 度 不 符 值 的 方差 降低 1,541.75 cms, 平均 降低 了 1.64 cm?， 对 应 于 1.28 cm 的 均 方 根 (RMS) 海 面 高 
C 沿 轨 海面 高 度 异 常 (Sea Level Anomalies, SLA) 也 是 估计 测 高 系统 整体 表现 的 一 个 重要 工具 , 党 轨 SLA 是 
N 通过 沿 轨 的 瞬时 海面 高 度 减 去 相应 的 平均 海面 得 到 , 本 文 所 用 到 的 平均 海面 由 CLS2011 平均 海面 模型 得 到 。 
v 沿 轨 SLA 的 方差 可 以 用 来 监测 测 高 系统 的 长 期 稳定 性 ， 其 值 越 小 表明 测 高 系统 的 整体 性 能 越 好 [is 202, JA E 
| 1 可 以 看 出 ， 和 Jason-e 的 GDR 中 的 SSB 校正 项 相 比 ， 在 其 它 条 件 都 相同 时 ，CR3D 可 以 将 沿 轨 SLA 的 方 

差 降 低 0.71~1.10 cm?， 平 均 降 低 了 0.92 cm?， 对 应 于 0.96 cm 的 均 方 根 (RMS) 海 面 高 度 。 

综 上 所 述 ， 和 Jason-2 中 的 SSB 校正 项 相 比 ，CRsD 可 以 大 幅 降低 海面 高 度 不 符 值 的 方差 和 沿 轨 SLA 的 

» 方差 ， 这 表明 CRaD 可 以 大 幅 提高 雷达 高 度 计 的 测 高 精度 
© 文献 [15,16] 的 3 维 SSB 模型 是 通过 海面 高 度 异 常数 据 得 到 的 ， 为 了 将 利用 本 文 方法 得 到 的 SSB 估计 值 
pm 和 利用 文献 [15,16] 的 3 4E SSB 模型 得 到 的 SSB 估计 值 进行 比较 ， 本 文 利用 2010 年 的 海面 高 度 异 常数 据 建 
Pe 立 了 一 个 3 维 SSB 模型 ， 并 使 用 该 模型 得 到 了 相应 的 SSB 估计 值 (SLA3D)。 图 4 给 出 了 相对 海面 高 度 不 符 

值 的 方差 ， 它 是 通过 分 别 利 ] CR3D 和 SLAsD 进行 SSB 校正 后 得 到 的 海面 高 度 不 符 值 的 方差 减 去 利用 

JA2GDR 进行 SSB 校正 后 的 海面 高 度 不 符 值 的 方差 得 到 。 图 5 给 出 了 相对 沿 轨 SLA 方差， 它 是 通过 分 别 利 
用 CR3D fll SLAsD 进行 SSB 校正 后 得 到 的 沿 轨 SLA 方差 减 去 利用 JA2GDR 进行 SSB 校正 后 的 沿 轨 SLA 77 
差 得 到 。 图 4 和 图 5 的 结果 都 是 利用 2011 年 (Cycle 93 - Cycle 127) 的 数据 通过 每 个 Cycle 单独 计算 得 到 。 从 
图 4 和 图 5 可 以 看 出 ， 和 JA2GDR 相 比 ，CRsD 和 SLAsD 都 可 以 降低 交叉 点 海面 高 度 不 符 值 的 方差 和 沿 轨 
SLA 的 方差 。 而 和 SLAsD 相 比 ， 利 用 CRsD 得 到 的 交叉 点 海面 高 度 不 符 值 的 方差 和 沿 轨 SLA 的 方差 都 更 
小 ， 这 也 表明 CRsD 比 SLA3D 有 更 好 的 性 能 
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6 结论 
海 况 偏差 由 于 其 较 大 的 不 确定 度 已 经 成 为 雷达 高 度 计 测量 海面 高 度 最 大 的 误差 来 源 。 传 统 的 海 况 偏差 


校正 采用 以 高 度 计 测 得 的 风速 和 有 效 波 高 为 输入 参量 的 2 维 经 验 模型 ， 本 文 从 非 线性 海浪 理论 的 角度 简 
高 差 的 大 致 关 系 ， 并 联合 ECMWF 再 分 析 数 据 集 提供 的 平均 波 周 期 以 及 高 度 计 测 得 
风速 和 有 效 波 高 ， 提 出 了 一 种 利用 交叉 点 数据 和 三 维 非 参数 模型 估计 海 况 偏差 的 方法 。 利 用 Jason-2 卫星 雷 
达 高 度 计 2009~2011 三 年 的 数据 对 得 到 的 SSB 估计 值 进行 评估 ， 


分 析 了 波 周 期 和 海 况 1 


图 5 相对 方差 (CR3D 和 SLA3D 得 到 的 SLA 的 方差 - JA2GDR 得 到 的 SLA 的 方差 ) 


Cycle 
4 ”相对 方差 “CRs3D 和 SLASD 得 到 的 海面 高 度 不 符 值 的 方差 - JA2GDR 得 到 的 海面 高 度 不 符 值 的 方差 ) 


单 
的 


和 Jason-2 的 GDR 产品 中 的 SSB 校正 项 相 


比 ,利用 本 文 方法 所 得 到 的 SSB 估计 值 能 将 交叉 点 海面 高 度 不 符 值 的 方差 和 沿 轨 SLA 方差 分别 降低 1.64 cm? 


和 0.92 cm?， 分 别 对 应 于 1.28 cm 和 0.96 cm 的 均 方 根 (RMS) 海 
到 的 s HÈ SSB 模型 相 比 ， 利 用 本 文 方法 所 得 到 的 3 28 


四 高 度 。 而 和 现 有 的 利 ) 
ESSB 模型 也 具有 更 好 的 表现 。 由 于 ECMWF 的 再 分 析 


海面 高 度 异 常数 据 得 


数据 已 经 被 广泛 应 用 于 业务 化 雷达 高 度 计 二 级 数据 处 理 ， 所 以 本 文 提出 的 海 况 偏差 估计 方法 具有 较 强 的 实 
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